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Abstract 

A nuclear density functional can be used to find the binding energy and shell structure of nuclei 
and the energy gap in superconducting nuclear matter. In this paper, we study the possible 
application of a nuclear density functional theory to nuclear astrophysics. From energy density 
functional theory, we can deduce the interaction between nucleons to find a rough estimate of 
the charge radius of the specific nuclei. Compared to the Finite-Range Thomas Fermi model, we 
include three-body forces, which might be important at densities several times that of nuclear matter 
density. We also add the momentum dependent interaction to take into account the effective mass 
of the nucleons. We study matter in the neutron star crust using the Wigner-Seitz cell method. By 
constructing the mass-radius relation of neutron stars and investigating lepton-rich nuclear matter 
in proto-neutron stars, we find that the density functional can be used to construct an equation of 
state of hot dense matter. 
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I. INTRODUCTION 



Understanding the equation of state (EOS) of dense matter is important to understand 
heavy-ion colhsion, supernova explosions, and neutron stars. Neutron stars are beheved to 
be composed of an outer crust, an inner crust, an outer core, and an inner core. The outer 
crust (10^ < p< 10" g cm ) has a lattice of neutron-rich nuclei in a gas of free electrons. 
In the inner crust (10^^ < P < 10^^ g cm~^), there are neutron-rich nuclei in a gas of free 
electrons and neutrons. The density of the outer core extends to ~ 5 x 10^'^ g cm~^ and is a 
homogeneous liquid mainly composed of neutrons, electrons, protons and muons. Not much 
is not known about the inner core (5 ~ lOpo) where the saturation density is po = 0.16fm^^), 
however, we believe that there are hyperons in the hadronic phase and deconfined quark 
matter 



24( 1 ■ A proper EOS should be able to explain nuclear properties in all density ranges. 
We can get basic information about nuclear matter from semiempirical mass formulas. One 
such fomula, known as the Liquid Drop Model (LDM), determines the binding energy(i?_) of 
nuclei in terms of a nuclear matter contribution and various corrections for finite nuclei|5|: 

(N — 

B{A) = a™iA + asurf^'/^ + a^oui^'^^'/^ + a^^m- — , (1) 

where a^i ^ -16MeV, Osurf ~ 20MeV, Ocoui ~ 0.751MeV, a^y.^ ~ 21.4MeV, is the number 
of neutrons, Z is the number of protons, and A is the total number of nucleons in a nucleus jsl. 
First is the volume term, which is the binding energy of infinite nuclear matter. The sec- 
ond term indicates the reduction of binding due to the nucleons on the surface. The third 
term represents the Coulomb energy, which is assumed to be that of uniformly-distributed 
charge. The last term is the symmetry energy, representing the decrease in binding energy 
for unequal numbers of protons and neutrons. 



Baym, Bethe, and Pethick (EBP) 26| made an EOS in the range of densities from 
4.3 X lO^^g/cm'^, where neutrons begin to drip out of the neuclei, to 5 x lO^^g/cm'^ (2po) us- 
ing a compressible LDM designed to take into account three important features: (i) the free 
neutron gas due to neutron drip; (ii) the nuclear surface energy reduced by the neutron gas ; 

T 1 

and (iii) the effect of the nuclear lattice Coulomb energy. Baym, Pethick, and Sutherland 27| 
extended the previous model in the regime from lO'^g/cm^ to neutron drip and applied it to 
neutron stars and white dwarfs. 



Oyamatsu|25|] studied nuclear shapes and lattice types (spherical, cylindrical, slab, cylindricl 



hole and spherical hole nuclei) at T = MeV with parametrized neutron and proton dis- 
tributions, and performed the Thomas-Fermi calculations with four different energy-density 
functionals. Negele and Vautherin 28| used a different method (Hartree-Fock approximation) 
to study the equation of state at T = MeV. They were able to get a density- dependent 
Hamiltonian from two-body interactions and found results similar to BBP. 
There were several efforts to make thermodynamic tables that can be used for supernova 
simulations. One of them is the Lattimer-Swesty EOS (LS)j3|. Their equation of state is an 
extension of the previous work by Lattimer et al.(LLPR) The LLPR EOS is based on 
the LDM. They included the effect of temperature on nuclei, the increase of surface energy 
as temperature increases, the effect of external nucleons, and the effect of nuclear excited 
states. The LS EOS also took into account nuclear deformation and the phase transitions 
from nuclei to uniform nuclear matter at subnuclear densities. Another table was given 
by Shen et al. 30|, in which they used the field-theoretical model with the Thomas- Fermi 
method. They constructed an EOS of nuclear matter in a wide range of the baryon mass 
density {ps = 1.25 x 10^ — 2.5 x 10^^ g/cm^), temperature (T = — 100 MeV), and proton 
fraction {Yp = — 0.56). 

Recently, Shen et al.[31] used a density- dependent relativistic mean-field theory to construct 
a nuclear EOS. For high- and intermediate-density nuclear matter, they employed relativisitic 
mean field calculations and used the virial expansion to study low-density nuclear matter. 
The table has the range of density from pB = 10^^ — 1.6fm~^, proton fraction Yp = — 0.56, 
and temperature T = 0.16 — 15.8MeV for high and low density nuclear matter. 
In this paper we use a simple density functional model to decribe both high and low nuclear 
density. To account for the short range of nuclear forces, we use a Gaussian form for the 
interaction. To find the minimum energy of a nuclear system, we use the Lagrange multiplier 
method; that is, the chemical potentials of neutrons and protons are constant in the cells of 
nuclear systems. Using this method, we can find the properties of single nuclei and heavy 
nuclei with a neutron gas in the neutron star crust. The MIT bag model is used to see 
the phase transition at high baryon density {ps > Po)- The mass-radius relation and the 
moment of inertia of cold neutron stars are calculated using the nuclear density functional. 
Proto-neutron star matter with neutrinos is also investigated for a given entropy per baryon. 
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II. NUCLEAR DENSITY FUNCTIONAL THEORY 



A. Energy density functional 



The energy of the nuclear matter can be given by 



E = Tkin + EpR + EzR + Ec + Es,L 



(2) 



where T^j^, EpR, Ezr, Ec, and Es^l are the kinetic energy, nuclear finite-range interac- 
tion, zero-range interaction. Coulomb interaction, and spin-orbit coupling respectively. The 
kinetic energy contribution from nucleons is simply obtained by 

^2 



(3) 



where t is the type of nucleons. 

In this Gaussian nuclear density functional (GNDF) theory, the number and kinetic densities 
are 



where ft is the Fermi-Dirac density function. 



1 



(5) 



'* 1 + e(^'-/*t)/'^ ■ 

For the finite-range term, we use a Gaussian phenomenological model for the nuclear poten- 



tial. 



+ 



(fri(fr2e ''12 /''o 



ViLpt{ri)pt{r2) + Viupt{ri)pt'{r2) 
V2Lpl^'{n)pl+^{f2) + V2upl^'{n)pmr2) 



(6) 



12 7^3/2^3 / 

+ J2 / d^ridh2e-'''^^''° J d^ptid^pt2ftift2V3Lpl2 + J d^ptid^Pt'2ftift'2V3upl2 

where pi2 = |pi -P2I, ?^i2 = |n - ^2!, tq is the length of interaction, and Vil, Vw, V2L, V2U, 
VsL, and V^u are interaction parameters to be determined. The last term is added to explain 
the effective mass of nucleons in dense matter. 

The zero-range term in the nuclear force can be regarded as the energy contribution from 
three-body nuclear forces. The three-body force is quite important if the baryon density 



increases beyond two or three times the saturation density. One possible form of the three- 
body force isjsl 

EzR = \h j d^rpn{r)pp{r)p{r) , (7) 

where is the interaction strength for three-body force, p„. (Pp) is neutron (proton) density, 
and p is total density. 

The energy functional for the Coulomb interaction has an exchange term which is absent in 
classical physics. 



Ec = E^J + E^f^ 



'2 J J ri2 
In bulk nuclear matter, the spin-orbit and Coulomb interactions constitute a small portion 
of the total energy, so we neglect these two terms in bulk nuclear matter. Then the bulk 
density functional would be 

Eb = Tkin + EpR + EzR 



where £b is the energy density for bulk matter. Using eq. (|3]), ([6]), and ([7j), we find 



(9) 
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B 



2m 2m 

+ ViL{pn{pn) + Pp{pp)) + Viu{pn{pp) + Pp{pn)) 

+ V2l{pI;''{pI;'') + pI^'{pI''')) + v^2r/(pi+^(pi+0 + pi+^(Pn+0) (lo) 

+ V3L(p„(r„) + Tn{pn) + Pp(l"p) + ^^(Pp)) + V3C7 (p„(rp) + r„(pp) + pp{T„) + Tp{pn)) 
1 

+ -^hPPnPp 

where we defined the Gaussian- type integral using '(...)': 

(^(^1)) = I d'r2e-^-^/Mr2) • (11) 

^ 

B. Effective mass, potential, and thermodynamic properties 

The effective mass of nucleons at the nuclear saturation density is about 0.7mB{mB = 
938)MeV. Some nuclear density functionals use the effective mass m* = mB] we, however. 
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introduced the momentum-dependent interaction, which describes the effective mass of nu- 
cleons. The functional derivative SEb gives us the effective masses and potentials in the 
nuclear density functional, 

6Eb = J d'r (vjpn + ^jTn + V,6p, + . (12) 

Now we get the effective mass for neutrons and protons. 



l + 4m{V,L{pt) + V,u{pt'))/n^ ^ ' 



and the potentials. 



ViLiPt) + VwiPt') + (1 + e)pl{V2L{pr) + V2u{pln) + VsLin) + Vsuirt') 



1 , 

+ -t3{2pt + Pt')Pt' ■ 



(14) 



where t' is a different type of nucleon from t nucleon. 

The thermodynamic properties are extremely important for describing the properties of hot, 
dense matter. The degeneracy parameter in the Fermi-Dirac distribution function is the key 

to the thermodynamic properties. Using Fermi-intcgrals, we get the baryon number density 
and kinetic density in terms of the degeneracy parameter (pt = [pt — Vt)/T, 

Landau's quasi-particle formula gives us the entropy density St, which tells us how to find 
the pressure in this density functional, 

-5. = / d'p[ft^nU + (1 - ft) ln(l - /,)] = ^^r, - ^^p, . (16) 
Prom the thermodynamic identity and entropy density given above, we can get the pressure: 

P = PnPn + PpPp + TSn + TSp - £ 

= Egn..,.)-. 

= + + yiL{pn{pn) + Pp{pp)) + Viu{pn{pp) + Pp{pn)) 

+ v,,{i + 26)(p^+^(p^j+^) + pj+^(p;+^)) + v,u{i + 26)(pi+^(p;+o + p'ripl^')) 
77 77 

+ V3L{pn{rn) + -^rn{pn) + Pp{rp) + -^Tpipp)) + V3u{pn{Tp) + :^Tn{pp) + pp{Tn) + ^^^(pn)) 



1 

+ ^hPPnPp , 

(17) 



where in zero-temperature, non-uniform matter, the chemical potential of the protons and 
neutrons are given by. 



2m; 

_^ 
2m 

+ 2 



(18) 



VlL{pt) + ViuiPt') + (1 + e)Pt(V^2L(p^^) + V^2t/(PD) + + ^3C/(Tt') 



1 



+ -t3{2pt + pt')pt' ■ 



III. PARAMETERS FOR THE GAUSSIAN NUCLEAR DENSITY FUNCTIONAL 



Every nuclear model should reproduce five nuclear matter properties: hinging energy, 
pressure, nuclear incompressibility, symmetry energy and efi^ective mass, m*. We use the 
saturation properties of nuclear matter to determine the parameters of the density functional. 
For zero-temperature, uniform nuclear matter, we have the energy density as a function of 
u = p/po and x = pp/p, 

= h'/^'/' [(1 - xf' + x'/'] + [v,,ix' + (1 - x)') + 2v,ux{l - x)] 
J-oPo t) 

+ 2'+'^u'+'^ [v2l{x'+'' + (1 - xf+'') + 2V2UX'+'{1 - xf^^] (19) 
+ 2'^\'/' [vsl{x'^' + (1 - xf') + vsu{x{l - xf/' + x'^'il - x))] + y.u'xil - x) , 

where we define the parameters 



V3L,U 



(20) 



5Tn 



Now we assume that the momentum-dependent interaction is blind to the type of nucleon, so 
VsL — Vsui^s — VsL + Vsu). The binding enery of symmetric nuclear matter {u — 1, x — 1/2) 
is then given by 



Po 



(21) 
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where Bq = 16MeV is the binding energy per baryon at the nuclear saturation density. 
The pressure at the saturation density vanishes, which mean the enery per baryon has its 
minimum at the saturation density, 



Po = PqTq 



^ + ^ + (1 + 2e)(f2L + V2U) + ^V3 + -tg 



0. 



(22) 



The incompressibihty parameter at the saturation density is given by 



dp 



P=PO 

265MeV . 



To [6 + 9iviL + viu) + 9(1 + 2e)(2 + 2e){v2L + V2u) + 20v3 + '^t'.^] 

o 



(23) 



The symmetry energy in nuclear matter is defined as 



8 dx^ 



p=p0,x=l/2 

I + ^^^^^ + (1 + + 2e)v2L - V2U) + - 



= 28MeV . 

Another parameter, which is related to symmetry energy, is given by 
^ 3pod^{S/p) 

p=po,x=l/2 



(24) 



8 dpdx"^ 



2 3 25 3 

- + -{viL - vw) + 3(1 + 2e)(l + e)((l + 2e)v2L - V2u) + —V3 - ^t', 



(25) 



= ^0 

= 54MeV . 

We choose the effective mass at the saturation density as 0.78m6 and use this number in the 



eq. ■ 



m 



m 



0.78mfe. 



(26) 



1 + 2mpoV^3/^' 

Thus we can easily recover ^3 from eq. fl20|) . From eq. fl2T]) . f l22|) . and f l23|) . we can have 
vi = viL + V2L, V2 = V2L + V2U and tg. 

5K0/T0 + 5^73(1 - 3e) - 72e - 905o/To(l + 2e) - 12 



V2 



Abe 

12 + 905o/To - 5Ko/To - bv^ 



90e(l - 2e) 



(27) 



to = —8vi — 16v2 — 8vs 



165o 48 



Tn 



5 
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TABLE I: Interaction parameters when e = 1/6, = 265 MeV, Sy = 28MeV, L = 54MeV. 



VlL 


viu 


V2L 


V2U 


V3L,U 


*'3 


ro (fm) 


-1.766 


-3.472 


0.410 


0.931 


0.169 


1.177 


1.205 



Then we can manipulate eq. ( |24|) . and (!25|) to get vn and V2l 

"5(1 -3e) 



f2L 



1 1 
-vi H 

2 2e 



27 



-t;3 



1 



-U2 



2(1 + e)'^^ 4e(l + e)2 
and we can have vm = f i — f il and V2u = f 2 — 'y2L- 



2e - 1 , (1 + 2e)S^ L 1 + 6e 
16 To 3Th 9~ 



5 ^ 1 

27^^ ~16^7^~37i^~9 



(2^ 



A. Determination of 1 + e power 



We added in eq (|6]) the auxiUary density interaction with the 1 + e power. We might 
regard 1 + e as the many-body effect — for example, a three-body force if e > |. It is known, 
however, that interactions among more than three-bodies are unimportant in dense matter. 
Thus we might restrict the e to be less than |. As the e changes, the ^3 parameter changes 
sign which means the three-body force can be attractive or repulsive. In the general Skyrme 
model with the three-body force, the parameter is positive. We choose e = 1/6 so that 
the interaction has the form of pJx Pt^ ■ ■'■^ zero-temperature, uniform matter, we have m^/^ 
terms in the energy density. From eq. f|T9|) . the energy density has u^/^, m^, u^/^, u®/^, and 

terms if we have e = 1/6 so we can use a statistical approach in uniform matter. Fig 1. 
shows the energy per baryon from GNDF and APR HI] EOS. We can see that as the density 
increases, the pressure from the two models agrees very well. 



B. The effective range of the nuclear force : tq 

In the Gaussian-interaction model, we can see the effective range of the force is given by 
ro, which is approximately ~ 1 to 2 fm. We don't have an analytic form of tq, so we need 
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-20 E 



0.00 0.05 0.10 0.15 0.20 0.25 



FIG. 1: The solid line represents the enery per baryon (uniform matter) using GNDF. The upper 
(lower) curve represents the energy per baryon of pure neutron matter (symmetric nuclear matter). 



The energy per baryon (dotted line) from the APR 



111 ] EOS was added for comparison. 



to rely on the numerical solution of the surface tension of semi-infinite nuclear matter: 

/oo /'OO 
\E - TSn - TSp - i2nPn " PpPp + Po]dz = - / [p{z) - po] dz , (29) 
oo J — oo 

where po is the pressure at 2; = — oo or 2; = +cxd. In one-dimensional, semi-infinite nuclear 
matter, we assume that the nuclear density depends only on the z-axis; the Gaussian integral 
then becomes 

—-3 / d^ru{r) = \ dzu{z) . (30) 

Experimental values of surface tension and surface thickness are u = 1.250MeV fm~^ and 
tgo-io = 2.3fm. Fig.2 shows the numerical calculation, which says that tq = 1.205fm from the 
surface tension and tq = 1.149 from tgo-io thickness. There is a 5% discrepancy between the 
two results. Table 1 shows the interaction parameters which we use in this paper when K = 
265MeV, Sy = 28MeV, L = 54MeV, and e = 1/6. The simple density dependent interactions 
(■^iL.c/) are attractive, on the other hand, the auxiliary density dependent interactions {v2l,u), 
momentum dependent ineractions {v3l,u) three-body force(t3) are repulsive in our model. 
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0.6 b d 3.5 r 




FIG. 2: The left figure shows the quantity po — pi^) at the semi-infinite nuclear surface when 
ro = 1.205fm and Yg = 0.5. The surface tension from this configuration is (x> = 1.250MeV fm~^. 
The right figure shows the surface tension (solid line) and tgo-io thickness (dashed line) as a function 
of tq. When tq = 1.205, tgo-io = 2.412fm. The surface tension and tgo-io thickness are both linear 
functions of rg. 



IV. NUCLEAR MATTER AND NUCLEI 



A. Specific heat 

The specific heat of uniform nuclear matter can be obtained by 

(31) 



dE 



For a non-interacting Fermion gas, the specific heat increases linearly with temperature. 
When the temperature is low enough, we expect that the specific heat of the nuclear matter 
tends to behave like a free Fermion gas. The specific heat formula for degenerate gas is given 
by0 

Cv = \m*kFklT, (32) 

where m* is effective mass of a nucleon, kp is the Fermi momentum, and fc^ is Boltzman 
constant. 

However, as the temperature increases, the non-linear behavior of the specific heat comes 
out so the degenerate gas formula is no longer valid, since the nucleons deep inside the Fermi 
surface are excited. 
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1.5 




T (MeV) 

FIG. 3: This figure shows the specific heat per nucleon of uniform matter for different densities. If 
the temperature is low enough, the specific heat behaves linearly with temperature. 



To calculate the _specific heat of uniform nuclear matter, we use the Johns, Ellis, and Lattimer 
(JEL) method |2l|, which enables us to get the pressure, energy density and entropy density 
for a given degeneracy parameter. Fig. 3 shows the specific heat per nucleon of uniform 
nuclear matter. It shows the linear relation between the specfic heat and temperature at low 
temperation as in eq. 

A detailed calculation of the specific heat at sub nuclear density in the neutron star needs to 
take into account the beta-equilibrium condition and heavy nuclei with a neutron gas. The 
specific heat plays an important role in the cooling process of neutron stars. In the neutron 
star crust, there are heavy nuclei and a free-neutron gas. The effective masses of protons 
and neutrons are different from the center of the heavy nuclei and dilute neutron gas, so we 
can't use eq. (15^ . In this case the specific heat at the neutron star crust can be calculated 
numerically by changing temperature and comparing the total energy change. 
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FIG. 4: This figure shows the basic properties of the closed-shell nuclei which can be obtained from 
the GNDF. The solid (dotted) line indicates the neutron (proton) density as a function of radius. 
As the number of nucelons in the nuclei increase, the neutron skin thickness(i?„ — Rp) increases. 



B. Nuclei at T = MeV 



We can use the GNDF theory and the Lagrange multiplier method to find the radius and 
binding energy per nucleon for a single nucleus using the Winger-Seitz cell method. In the 
Lagrange multiplier method, the chemical potentials of protons and neutrons are constant 
in the Wigner-Seitz cell to minimize the total free energy. Fig. 4 shows the radius and 
binding energy of the closed-shell nuclei using this method. These results agree well with 



experiment 



131]. ^°Ca has a larger charge radius than neutron radius because of the Coulomb 



repulsion between protons. The solid (dotted) line denotes neutron (proton) density. As 
the atomic number increases, the central density of neutrons increases; on the other hand 
the central density of protons decreases. The difference between charge and neutron radii 
increases and the neutron skin becomes thicker as the atomic number increases. In ^'^'^Pb 
nuclei, the central density of protons is lower than the proton density of the outer part of 
nuclei (r = 4 — 5 fm) because of proton Coulomb repulsion. 

Table II shows the proton and neutron radii and binding energy per baryon of closed shell 
nuclei from various nuclear models. The calculation from GNDF theory agrees well with 
experimental results. 
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TABLE II: Comparsion of the results from Steiner (Potential &: Field Theoretical) et al. 13j] , FRTF 
l|l|], FRTF II [2] and the GNDF 



Nucleus 


Property 


Experiment Potential 


FT 


FRTF I FRTF II GNDF 


208 p]^ 


rch (fm) 


5.50 


5.41 


5.41 


5.38 


5.45 


5.44 




BE/A(MeV) 


7.87 


7.87 


7.77 


8.01 


8.17 


8.02 




5R(fm) 


0.12 ± 0.05 


0.19 


0.20 


0.15 


0.13 


0.14 






0.20 ± 0.04 












90Zr 


rch (fm) 


4.27 


4.18 


4.17 


4.10 


4.15 


4.17 




BE/A(MeV) 


8.71 


8.88 


8.65 


8.77 


9.00 


8.72 




5R(fm) 


0.09 ± 0.07 


0.075 


0.093 


0.064 


0.054 


0.057 


40Ca 


rch (fm) 


3.48 


3.40 


3.34 


3.22 


3.26 


3.31 




BE/A(MeV) 


8.45 


8.89 


8.61 


8.47 


8.77 


8.33 




5R(fm) 


-0.06 ± 0.05 


-0.044 


-0.046 


-0.036 


-0.039 


-0.042 



-0.05 ± 0.04 



C. Heavy nuclei in the neutron star crust 

In the neutron star crust, heavy nuclei are formed with a free- neutron gas. These heavy 
nuclei are suspected to form a BCC (body centered cubic) structure. In the static equilibrium 
state, we calculate the density profile of heavy nuclei with a neutron gas using the Wigner- 
Seitz Cell method. The plot on the left side of Fig. 5 shows the proton (dotted line) 
and neutron density (solid line) profiles from the center (r=0 fm) of heavy nuclei when 
p = O.Olfm"^. There are dripped neutrons outside of the heavy nuclei. The cell size (Rc), 
which is a rough estimate of the distance between neighboring heavy nuclei, is determined 
by nuclear density and beta equilibrium conditions (/i„ = /ip + /ie). There is a wave function 
overlap at the boundary of Wigner-Setiz cell. The actual distance between heavy nuclei is 
(87r/3)^/'^Rc. The right side of Fig. 5 shows the binding energy per baryon as a function of 
Wigner-Seitz cell size. As the density decreases, the cell size increases and the energy per 
baryon converges to -8.0 MeV. 

Table III shows the thermodynamic properties and physical dimensions of nuclei in the 
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0.3 
0.2 
0.1 
0.0 





R, = 29.155fm ; 




N =1003.7 J 




N„„c = '^7.7 




\ Z = 34.3 J 




\ n = 0.010/fm^ ; 




\ E/A = 2.21 8MeS 




\ = 2.452MeV ; 




\ R„ = 6.200fm J 




\ R„ = 4.960fm ; 











(fm) 



10 12 




200 400 600 800 1000 1200 1400 
(fm) 



FIG. 5: In a neutron star, heavy nuclei exist. The left figure shows the density profile of proton 
(dotted line) and neutron (solid line), r = fm means the center of heavy nuclei. Outside the heavy 
nuclei, there are dripped neutrons. As density decreases, the cell size increases and the energy per 
baryon converges to -8.0 MeV. 



TABLE III: Nuclear properties in the neutron star crust 



p (fm^) 




p (MeV/fm^) 


e (MeV/fm^) 


Nnuc 


Z 


Rc (fm 


5.623x10- 


-2 


0.181 




53.06 




271.6 


89.51 


25.12 


5.012x10" 


-2 


0.147 




47.27 




218.0 


32.18 


18.70 


3.981x10" 


-2 


9.332x10" 


-2 


37.52 




137.8 


21.69 


18.06 


2.512x10" 


-2 


3.965x10- 


-2 


23.65 




145.3 


27.93 


22.77 


1.585x10" 


-2 


1.976x10" 


-2 


14.91 




150.7 


32.07 


26.36 


1x10-2 




1.081x10" 


-2 


9.405 




147.7 


34.29 


29.15 


1x10-3 




7.988x10" 


-4 


9.383x10- 


-1 


116.3 


38.50 


40.71 


1x10-^ 




6.945x10" 


-5 


9.352x10" 


-2 


79.23 


38.31 


65.47 


1x10-^ 




1.480x10" 


-6 


9.326x10" 


-3 


58.02 


36.98 


131.4 


1x10-6 




2.362x10" 


-9 


9.311x10" 


-4 


47.90 


35.14 


270.6 


1x10-^ 




1.714x10" 


-9 


9.303x10" 


-5 


43.37 


34.10 


569.7 



neutron star crust. Nnuc and Z are number of neutrons and protons of heavy nuclei in the 
Wigncr-Seitz Cell. The atomic number of heavy nuclei remains Z ~ 35 for a large range of 
densities before the phase transition to uniform matter. This means that the proton fraction 
decreases as the density increases. For a narrow range of densities, the atomic number 
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suddenly increases and the heavy nuclei merge with free neutrons to make uniform nuclear 
matter. 
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FIG. 6: Effective mass of nucleons in the Wigner-Seitz cell as a function of radial distance from 
the center of heavy nuclei. Since the nuclear interaction is weak at the boundary of the cell, the 
effective mass of nucleon and the pure mass of a nucleon become equal. 

Fig. 6 shows the effective mass of nucleons in the Wigner-Seitz cell. The effective mass of 
nucleons in the Wigner-Seitz cell is given by eq. (fT^ and eq. f EUj) . 



Since we assume the momentum interaction is blind with respect to isospin, the effective 
mass is identical for different isospin nucleons. The effective mass to pure mass ratio of 
nucleons is 0.78 at the center of the heavy nuclei and becomes 1 at the outer region of the 
Wigner-Seitz cells since the density of nuclear matter is low at the outside of the heavy 
nuclei, the interaction energy of nuclear matter is weak. 

V. PHASE TRANSITION 

In the neutron star, we can see two types of phase transitions: one is the phase transition 
from nuclei with a neutron gas to uniform matter, and the other is the phase transition from 



m 



(33) 
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uniform nuclear matter to quark matter. During the first phase transition, we can see the 
nuclear pasta phase. That is, spherical nuclei become ellipsoidal, then cylindrical, and finally 
slab phase before nuclear matter becomes uniform matter. However, the energy difference 
is quite small, so that the effects on the large scale physics are neglible. On the other hand, 
the second phase transition is quite dramatic. The energy and pressure change significantly 
from nuclear matter to quark matter. 



A. Uniform matter 



To check the phase transition points from heavy nuclei with a neutron gas to uniform 
nuclear matter, we can simply compare the energy per baryon of uniform nuclear matter 
with the energy per baryon of nuclei with a neutron gas since the nuclear matter exists in 
the lowest energy states. The energy per baryon in uniform matter can be easily obtained 
by changing the '(...)' integrals to non integral form from eq. ( lT2l) since the Gaussian inte- 
grations in uniform matter become unity. Typically there is a phase transition around 0.5po- 
We know that in the outer crust of the neutron star the nuclei has a BCC structure. If we 
assume that the pasta phase exists in the low-density region, we may use the density pertur- 
bation to see the phase transition from nuclei with a neutron gas to uniform nuclear matter. 
We use the wave number pertubation to see the energy exchange, which has contributions 
from the volume effects, gradient effects, and Coulomb energy can be approximated [K 



v{q)^v, + Pq^+ , (34) 

q -\- fhrpp 

where q is the sinusoidal variation of the wave number in the spatially periodic density 
perturbation. 

The volume term is given by 

° dpp (dfln/dpn) 

The energy exchange from the gradient has the form 

/3 = D,, + 2Dr.,^ + DU\ ^ = -P^ (36) 

where the coefficients of the gradient terms are given by Dpp = Dnp = Dnn = 132MeV ■ 
fm^ jl?!]. The kxp in the Coulomb interaction represents the inverse Thomas- Fermi screening 
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length of the electrons. When we see the change in the sign of v, the uniform matter phase 
is more stable than the perodic structure of the nuclei. The v has a minimum at 



Vmin = Vo + 2(47^6^(3)'/^ - f3k 



TF 5 



(37) 



when q. 



2 

min 



TF- 



Another way to see the phase transition is to use the thermodynamic instability. The ther- 
modynamic stability condition can be described using the inequalities jl^ jl?! . 



dP 

dv 

dqc 



>0, 



>0. 



(38) 



where P = Pb + Pe is the total pressure from electrons and baryons and /i = fin — fJ'p is 
the difference between the neutron and proton chemical potentials, which is the electron 
chemical potential in beta-stable matter. Qc is defined clS — Xp Pg / p. Mathematically, 
the inequalities in eq. ( 138|) show that the energy per baryon is convex. Eq. (!38|) can be 



verified to be 



zap 

\ dv 



p 



2p 



dE{p,Xp) ^d'^E{p,Xp) 



dp 
dqc 



dp 

d'^E{p,Xp) 



dx"^ 



+ 



dp^ 



d'^E{p,X p) ^\ ^ I d'^E{p, Xp 
dpdxp 



-P 



dx"^ 



> 0, 



(39) 



pI 



Ti'^h^p 



> 0. 



The second of eq. fl39l) always holds, so the first will determine the phase transition in the 
neutron star crust. In Xu et al.|17|], they use a simple equation to determine the instability 
using the thermodynamic relation. 



2dEd^E d^Ed^E / d'^E 
p dp dXp dp"^ dXp \dpdxp 



dpn dpp [ dpn 



dpn dpp 



dpp 



(40) 



Eq. ( HOj) is equilivalent to the volume part of the thermodynamic perturbation eq. (!35|) 
method except that there is a dpn/dpn difference. Comparing the two methods (pertubation 
and thermodynamic instability) shows the effects of the gradient and Coulomb terms in the 
pertubation method on the transition densities. Fig. 7 shows transtion densities using the 
perturbation method and thermodynamic instability. The perturbation method has a lower 
transition density (0.355po) than thermodynamic instibility method (0.406po)- 
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FIG. 7: We can see the transition density from nuclei with a neutron gas to nuclear matter. The 
solid line denotes the curve from the perturbation method. The dotted line corresponds to the 
thermodynamic instability method, equilivalent to vq in the perturbation method. The solid line 
has a transition density of 0.355po and the dotted line has a transition density of 0.406/jo- 

B. Quark matter 



In this paper, we don't consider the appearance of hyperons since it is not clear how the 
hyperons and nucleons interact. Thus we simply consider the phase transition from uniform 
matter to quark matter. We use the MIT bag model for the quark matter equation of state. 
At T=0 MeV, the pressure and energy density are given by [71] 



p 



I 2 2 4\l/2/ 2 ^ 2 4\ -'- 4 8 i 

^A^f - ) ' if^f- 2^fC ) - ^rrifC In 



(41) 



where the density for each quark flavor is given by 

_ (/i2 - m)fl^ 
7r2(nc)3 ■ 

For the pure-quark phase we use = rud = 0, rus = ISOMeV and B = lOOMeV fm"^. In 



(42) 
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FIG. 8: The left panel shows the energy density of nuclear matter (dotted line), quark matter 
(dashed) and mixed phase (solid). The right panel shows the pressure of nuclear matter (dotted 
lined), quark matter (dashed) and mixed phase (solid). 

the mixed phase of uniform nuclear matter and quark matter, we apply Gibb's conditions 
to minimize the free-energy density with two constraints, which are related to total number 
denstiy and charge neutrality. 



where x is the volume fraction of the uniform nuclear matter in the mixed phase and the 
subscript N (Q) represents nuclear (quark) matter. In the mixed phase, the total charge is 
globally neutral in contrast to pure nuclear matter and quark matter. Prom minimizing the 
free enegy, we have 



Pb = XPn + (1 - x)Pq 

<5 = xQiv + (i-x)gQ = o, 



(43) 



Pn^Pq 



flu + 2fXd 



(44) 



/ip — 2/i„ + /Id , 
then we have the energy density of the mixed phase 



e = X^N + (1 - x)eQ ■ 



(45) 
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FIG. 9: The left panel shows the mass-radius relation of the neutron stars (dotted line) from the 
GNDF and the hybrid stars (solid line). The right panel shows the number density profile of a 
hybrid star which have M= 1.44M0. At r ~ 5.0km, the phase transition to nuclear matter takes 
place, and protons and neutrons exit. When r=8.2km, the phase transition to nuclear matter is 
completed and there is no more quark matter. 

As in the case of uniform matter, we assume beta-stable matter so that the chemical poten- 
tials of the nuclear matter and quark matter have the relation 

(46) 

Fig. 8 shows the energy density and pressure as a function of number density. The dotted 
(dashed) line represents nuclear (quark) matter. The solid line denotes the mixed phase. The 
phase transition begins when the baryon density becomes 1.386po and all nucleons turn into 
quark matter when the baryon density becomes 5.236po- 

If there is a phase transition in the core of a cold neutron star, the mass and radius are 
quite different from the case of a pure-nuclcar-mattcr neutron star. The left panel of Fig. 
9 shows the mass-radius relation of the neutron stars (dotted line) and hybrid stars (solid 
line). The right panel shows the number density profiles of a quark matter and nuclear 
matter of 1.44M0 hyrid star. The maximum mass of a cold neutron star with mixed phase 
(hybrid star) is 1.441Mo and the central density of the neutron star is 9.585po- The mass 
and radius curve with the mixed phase indicates when the mixed phase happens. As the 
distance from the center of the hybrid star increases, the phase transition to nuclear matter 

21 



takes place so that neutrons and protons appear(r ~ 5.0km). Far from the center, the phase 
transition is completed (r = 8.2km); pure nuclear matter exists only for larger radii. The 
existence of quark matter in the hybrid star can be explained by angular momentum loss of 
the proto-neutron star. That is, fast-rotating neutron stars lose angular momentum because 
of magnetic dipole radiation; the central density of the neutron star increases due to the 
decrease in centripetal forec, then quark matter appears. Since quark matter has a lower 
energy density than pure nuclear matter, we might expect heating of the neutron star from 
latent heat from quark matter. 



VI. ASTROPHYSICAL APPLICATION 

A. Mass-radius relation of a cold neutron star 



We know that the radius of a neutron star is ~ 10km and the mass is ~ IAMq. In 
this system, the degeneracy pressure of the neutrons provides support against gravitational 
collapse. We can apply our model to calculate the mass and radius of neutron stars for 
a given central density. We use the Tolman-Oppenheimer-Volkov (TOV) equations which 
describe general relativistic hydrostatic equilibrium: 

dp G{M{r) + Anr^p/c^) (e + p) 

d^ ^ ^^(7^^2GM(07c^)c^ ^^^^ 
dM 

dr & 

Fig. 10 shows the mass-radius relation for a cold neutron star. In the GNDF model, the 
maximum mass of a cold neutron star is 2.163Mq, and the correspoding radius is 10.673km. 
The maximum mass from the GNDF model is in between the FRTF truncated model (FRTF 
I) and the modified model of the FRTF (FRTF II)[l||2|. 



B. Moment of inertia of the neutron star 



In the slow-motion approximation, the moment of inertia is given by|15| 

[''r\p + p)e'-^-'^/^oodr, 
3 Jo 



(4J 
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FIG. 10: Mass-radius relation for a cold neutron star. The mass of a cold neutron star from the 
Gaussian density functional model has a maximum mass of (2.163Mq) when the central density is 
6.74po. 



where A = — ln(l — 2m/r) and u are the metric coefficients and u is the rotational drag 
function. In terms of the function j = e'"^^^"^^'^ , the rotational drag statistics 



dr \ dr J dr 



(49) 



with the boundary conditions 



UJR=1- 
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dr\ °' 



(50) 



Therefore, the moment of inertia can be written as 



- / r uj — dr = - d\r 00 
3 Jo dr 6 Jo V dr 



4 .dj 



B^du 
6 dr 



(51) 



R 



We note that the second-order differential equation that u satisfies, eq. (149|) . can be instead 
written as a ffist order differential equation in terms of the function (p = dlnco/dlnr, 

d In J 



-^^ = -^(0 + 3)-(4 + ,, 
dr r dr 



where 



ciln j 

dr 



r — 2m 



(P + P) 



(52) 



(53) 
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with the boundary condition 0(0) = 0. The moment of inertia becomes 



o 



(54) 



using the boundary condition for u. This simphfies to 

I 



6 + 2(j)R 



(55) 



Lattimer and Schutz proposed an empirical approximation for the moment of inertia 22] , 
/ ~ (0.237 ±0.008)Mi?2 



Mkm /Mkm\ 
1 + 4.2^^^ + 90(T^^j 



(56) 



MqR \MqR. 

Fig. 11 shows the moment of inertia of a cold neutron star. The color band represents upper 
and lower boundaries of the emprical approximation. FRTF I, II and GNDF agree quite 
well with this empirical approximation. 
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FIG. 11: Moment of inertia of a cold neutron star. The color band represents the upper and lower 
boundaries of the emprical approximation eq (I56p . Three different models show different curves; 
however, they represent the emprical results quite well. 
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C. Matter in the proto-neutron star 



We expect that the core of a proto-neutron star is lepton rich (Yi ~ 0.4) and the entropy 



per baryon is s ~ 2 — 4 



23| . Since the Gaussian density functional model deals with non- 



relativistic nucleons, we use the relativistic leptonic formula to find the contribution from 
leptons to the total pressure, energy and entropy. Since the electrons and neutrinos are 
relativistic in proto-neutron star matter, we can get analytic solutions for leptons which 
are! 



r — 



/3 



(57) 



where t = 37i^{hc)^nYi/ gi and q 



TY — ^mfc'^. Here gi is the spin degeneracy; thus 



ge = 2, and g^ = 1 since the neutrino is only left handed. We assume that the neutrino is 
massless. 

Expressions for the pressure and entropy per baryon are 



Pi 



"247r2 \hc) 



6n{hcy 



1 + 



2rr2 ^ 2 4 
TT i m, C 



^ 2rr2 ^ 2 i- 

— TT I m, c 

15 2 ' 



(58) 



1 

2' 



To simplify the model of the core of a proto-neutron star, we assume that there are only 
protons, neutrons, electrons and electron neutrinos . In this case 

5^ ~t~ '^'jD ~t~ 5e ~1~ Sy 

(59) 

If the core of the proto neutron star is in equilibrium (to be more exact, it is in quasi-static 
equilibrium), then the EOS needs to meet eq. (158!) and eq. (1591) . On the other hand, in the 
cold, catalysed neutron star matter, we can simply say that there are only electrons. Thus 

Hn = jJ'e + fJ'p ■ (60) 
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FIG. 12: The left figure shows the total pressure when s=2.0, s=0.5, T=OMeV and T=15MeV. 
The total pressure from various conditions shows different aspects. The total pressure for a fixed 
entropy per baryon (s=2, s=0.5) with fixed lepton fraction (Yi=OA) shows identical results for 
lower densities, but they have different behavior as the density increases. The temperature for a 
fixed entropy per baryon can be obtained from beta-equilibrium. The right figure shows the lepton 
fraction for four different cases. The results from fixed s=2 and s=4 are not distinguishable. The 
lepton fraction for T=0 and 15MeV, we assume beta-equilibrium. 

If the interation of leptons with nuclear matter is weak, then the chemical potential and 
pressure of the leptons at zero temperature are given by[3|, 

-\ + -m^c*^ In f 
4/8 V mc 



2 



(61) 

In Fig. 12, the pressure is indistinguishable when s=2 and s=0.5 with fixed lepton fraction 
(Y; = 0.4) at lower densities. However, for higher densities, they show a different behavior. 
At low densities, the pressure at high T is greater than the ones from the fixed entropy 
per baryon with fixed lepton fraction, however, at high densities, the pressure contribution 
from the trapped neutrinos is much more important than the thermal effect. The right side 
of Fig. 12 shows the lepton fraction from electrons and neutrinos for a given entropy per 
baryon. The lepton fraction for each case is almost identical. For the cases T=0 and 15 
MeV, we assumed the beta-equilibrium condition. The electron fraction when T=15MeV 
increases from 0.06 to 0.25 between p = 0.3po < P = lOpo- Fig. 13 shows the various 
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FIG. 13: Entropy contribution from baryons and leptons for each fixed entropy per baryon. Since 
the neutron (electron) fraction is always greater than the proton (neutrino) fraction, the entropy 
contribution from neutrons (electrons) is greater than protons (neutrinos). 

entropy contribution for a given entropy per baryon and lepton fraction. Each graph shows 
the same trends. The entropy contributions from baryons is decreasing as baryon density 
decreases, and the electron's (neutron's) contribution to entropy is always greater than the 
neutrino's (proton's) contribution since the electron (neutron) fraction is greater than the 
neutrino (proton) fraction. 



VII. CONCLUSIONS 

In the GNDF model, the total energy consists of the kinetic energy, the finite-range 
effect, the zero-range effect, and the Coulomb energy. Using the Lagrange multipher 
method, we find the potential energy, pressure, chemical potential and thermodynamic 
properties. The interaction parameters were obtained from the properties of infinite nuclear 
matter. We can find the charge radius and binding energy per baryon of the closed-shell 
nuclei. The Wigner-Seitz cell size increases as the density decreases and the binding energy 
per baryon approaches -S.OMeV. The effective mass becomes 0.78m at the center of the 
heavy nuclei and becomes m outside of heavy nuclei. We were also able to find the pressure 
of uniform symmetric nuclear matter and neutron matter. For finite temperature, we can 
see the specific heat of nuclear matter follows the general trend of the free fermions. Thus 
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the GNDF model is a good nuclear matter model to study for both low and high nuclear 
densities. We can improve the current model if we have more exact experimental results 
and we add additional interaction terms to explain the experimental results. The phase 
transition was studied using the GNDF model. The density perturbation suggested that 
the phase transition from non-uniform nuclear matter to uniform matter takes place at 
densities less than 0.5po- When we take into account the phase transition from uniform 
nuclear matter to quark matter, we see there is a drastic change in the maximum mass 
of a neutron star, since the pressure and energy density of quark matter are significantly 
different from nuclear matter. The maximum mass of the hybrid star is less than 1.5Mq. 
This might rule out the coexistence of nuclear matter and quark matter in compact stars. 
The GNDF model can be used to study proto-neutron star matter combined with a leptonic 
environment. Using the GNDF theory, we may compare the maximum and minimum 
masses of proto-neutron stars. 

In a subsequent paper, we will study the nuclear pasta phase (spherical shell, cylinder, and 
slab geometry) using GNDF. 
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